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Abstract 

We check the recently proposed higher loop Bethe-ansatz for the sl{2) sector of A/" = 4 at 
two loops by a direct perturbative calculation using M = 2 superfields in super symmetric 
dimensional reduction. Our method can in principle address operators built from many 
elementary fields and of arbitrary twist. We work out in detail the spin three, twist three 
operator mixing problem at leading order in N and find agreement with the prediction 
based on integrability. 



1 Introduction 



The original version of the AdS5-CFT4 correspondence is a weak-strong duahty so that 
a direct comparison of string theory with gauge theory quantities is hardly possible in 
the perturbative regime. In a special limit of the geometry this difficulty may have been 
overcome by the BMN-proposal 1 . The exact leading orthogonalisation of the relevant 
set of operators has been connected to integrable spin chains, thus opening the possibility 
of analysing the large N mixing in = 4 through the Bethe-ansatz [^. The spin chain 
Hamiltonian can be identified with the dilatation operator of the planar CFT 151, a powerful 
tool that has successfully been used e.g. to compute one- and two-loop scaling dimensions 
[1] of the gauge theory equivalent of a class of string theory solitons with large angular 
momentum jSj. 

On the field theory side, dilatation operator and Bethe-ansatz are known for the full 
theory at one-loop jB|, while the dilatation operator has been worked out up to three loops 
for an sn(2|3) subsector and if the subalgebra is further restricted to su{2) also to four 
loops [SI and five loops [Hj. For the su{2) sector a Bethe-ansatz up to three loops was 
formulated in jJOl and there is even an all-loops conjecture based on a new long-range spin 
chain j^]. 

A step towards the quantisation of AdSsxSs string theory has been made in by 
taking into account the lowest term in a perturbation series in the curvature around the 
pp-wave background underlying the BMN construction. Data from here as well as from 
the soliton calculations has been incorporated into a Bethe-ansatz on the string side |12j . 
The ansatz deviates from the gauge theory conjecture starting at three loop order in the 
gauge coupling, where there are extra scattering terms, although these are expressed in 
terms of quantities known from the gauge theory [Sj. The three-loop disagreement of the 
Bethe-ansaetze accurately reproduces the gauge theory/string theory mismatches found up 
to now ^1^^. Despite of this, the integrable structures on the two sides of the conjecture 
are very similar: one may for example match the higher charges |13j at one and two loops. 
Furthermore, the spectrum of the classical string in the su{2) sector was derived from a set 
of integral Bethe-equations in the equivalent of the thermodynamic limit of the discrete 
Bethe equations in the gauge theory picture. The idea behind the string Bethe-ansatz was 
in fact to look for an appropriate discretisation of the continuum equations. The ansatz 
has the feature of "factorised scattering", i.e. that the Bethe equations involve a product 
of two-site interactions termed S'-matrices. 

The article shifts the attention entirely to the S-matrix. It is observed that the 
energy eigenvalues found in jJJ for an su(l|l) part of the sm(2|3) sector and in j^S] for 
the disjoint sl{2) sector are compatible with factorised scattering. A relation between the 
string theory 5-matrices in the three sectors is transferred to gauge theory, yielding a two- 
and three-loop Bethe-ansatz in the sl{2) sector. The associated spectrum is checked for the 
twist two operators, whose anomalous dimensions up to three loops were extracted in ^\ 
from an outstanding QCD calculation P^ . 

The purpose of the present paper is to test the two-loop anomalous dimensions (i.e. 
order g"^) given in ^Sj for the twist three, spin three operators in the sl{2) sector against a 
direct calculation of two-point functions based on AA = 2 superfields ^H] and regularisation 
by supersymmetric dimensional reduction (SSDR) [10].^ The result strongly vindicates the 

^Superspace calculations package ordinary Feynman graphs into supergraphs. In many situations this 
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iS-matrix conjecture of since at least in the dilatation operator picture there is the issue 
of "wrapping", i.e. how the general interaction is modified when there are only few sites 
in the spin chain. We emphasize that our method is applicable beyond twist two; barring 
for computing limitations we can handle operators built from many elementary fields and 
carrying higher spins. 

We encounter Feynman graphs with up to four loop integrations in momentum space. ^ 
The only four-loop diagram is of a type that we call "BPS-like" since the same supergraphs 
without extra partial derivatives on the outer legs occur in the two-point functions of 1/2 
BPS operators. For the two- and three-loop diagrams of this type we verify finiteness 
with and without extra partial derivatives using the cutting edge Mincer algorithm I 21j . 
By analogy, we assume that also the four-loop BPS-like supergraph does not contribute. 
On the other hand, all the singular integrals have maximally three loop integrations in 
momentum space so that they can be evaluated exactly. 

Note that the Mincer programme itself calculates massless two-point integrals with two 
or three loop integrations and arbitrary scalar numerators. We reconstruct integrals with 
open indices from their projections with the total ingoing momentum. 

Future work aims at the construction of the planar two-loop dilatation operator for 
the sl{2) sector, which would allow one to further analyse the recently discovered one-loop 
gauge/string theory discrepancies in the sl{2) sector (22] • 

2 The mixing problem 

We consider the leading N mixing of the operators 

{^1,^2,^3} = rr(Pf Z)(Z)f Z)(Z)f Z)) . (1) 

Here Z is a complex scalar field of the = 4 SYM theory with SU{N) gauge group and 
Dfj, = + ig is the Yang-Mills covariantised space-time derivative. The operators 
carry spin s = si + S2 + S3, i.e. the Lorentz indices are symmetrised and traceless, which 
can be made automatic by projecting all indices onto z = xi + 1x2 as indicated in 

The theory is regularised by SSDR (super symmetric dimensional reduction) 120] in x- 
space. This scheme is related to dimensional regularisation, where the scalar propagator is 

(Z(1)Z(2)) = 4(^x?2r, CO = □i(Z(l)Z(2)) = 5(xi2). (2) 

X12 

(The normalisation change involving non-integer powers of vr and ^(2) is hidden in the 
mass scale. We will not usually write //, since it can easily be reinstated.) On dimensional 
grounds, Feynman-graphs must come with the fractional powers 

9' : {xl,r, 9^ : {xl,r, 9' : {xl^f^ , (3) 

whose e-expansion yields logarithms. The conformal invar iance of A/" = 4 SYM implies that 
these logs exponentiate to scaling dimensions associated with composite operators. 

leads to a significant decrease of tfie number of diagrams. We prefer the M = 2 over the N = 1 formulation 
because the extra supersymmetry helps to suppress even more graphs. 

■^In the literature the 0{g^), 0{g'^) parts of the scaling dimension are usually called one-loop and two-loop 
anomalous dimension. Nevertheless, when the first two anomalous dimensions are computed from two-point 
functions the momentum space Feynman integrals have two loops at order and up to four at g'^. 
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Such conformal properties become visible only in the renormalised theory, i.e. after 
introduction of 2-factors and in the limit e — > 0. The set of operators has to be orthog- 
onalised, upon which the two-point function of a correctly renormalised primary of spin s 
has the form 

where the sets of indices {fii . . . fig} and {i^i . . . z/^} are both individually symmetrised and 
made traceless |23| . In the last equation we have used the inversion tensor 

Jfiuix) = rif,^ - 2^^^^ (5) 

to organise the terms in a simple manner. The scaling dimension A = Aq + (7'^7i + 5^72 + • • • 
splits into a classical part Aq and the "anomalous dimension" depending on g^. It can clearly 
be read off from any term in We find it convenient to drop all traces because this reduces 
the set of graphs, see below. If desired, the full functional form of any two-point function 
can be reconstructed: for primary operators we may use (jU while two-point functions of 
conformal descendants are obtained from it by differentiation. 

We are interested in the anomalous dimension up to order g^, which requires the leading 
and the next-to-leading term in the e-expansion, i.e. 

5° : 0(1), 0(e), g^ : O(-), 0(1), / : 0(1), O(-) . (6) 

e e 

In particular, at g^ we do not need to keep track of the finite part of the integrals. Cor- 
respondingly, we will not eliminate finite 0{g'^) cross terms between the operators. Terms 
of order e and higher yield contact terms when the regulator is taken off [251 • Although 
we are generally not interested in the contact part we cannot discard the 0(e) terms in 
the tree-level mixing, because the singular renormalisation factors relate them to the sub- 
leading terms at g'^ and g'^. The renormalisation/orthogonalisation proceeds along these 
lines: first, we need to choose a basis in which the 0(1) tree-level mixing and the simple 
logarithm terms at g^ are diagonalised. In such a basis each operator picks up a ^-factor 
of the standard form |2H] 

and 71,72 can already be calculated. Second, in order to fully orthogonalise we may have 
to subtract the other operators with mixing coefficients of the type 

B = g'ho + g'(^^+b2o)+0{g'), (8) 

(In this letter 620 remains undetermined because we truncate the e-expansion.) This form 
of the ^-factors may be called "minimal subtraction" since we strictly eliminate only sin- 
gularities and finite cross terms. 

At spin we have the operator 

O = {0,0,0}. (9) 
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At spin 1 we find its descendant 



Z),0 = {1,0,0}. (10) 

At spin 2 there are 

L>^0 = {2, 0,0} + 2 {1,1,0}, (11) 
Kg = {2,0,0}-(2 + e){l,l,0}. 

The operator O is a 1/2 BPS state. It is gratifying to see that its g^,g'^ two-point 
functions do indeed not contain singular or finite pieces. For xi ^ X2, e — > the loop 
contributions to {00) converge to zero and thus our calculation confirms protectedness. 
Let the descendants Df,0, s = 1,2,3 be defined by naive distribution of each derivative 
over the three sites in the operator. All three descendants are exact null-vectors of the 
and g^ mixing. 

At spin 2 there is a second operator which we must orthogonalise w.r.t. D^O. The 
definition (jllj) for Kq is not quite in the spirit of minimal subtraction because we have used 
a mixing coefficient involving a positive power of e. It has the advantage that the bare Kq is 
tree-orthogonal to D^O at leading and next-to- leading order in e, whereas this is guaranteed 
at g^,g^- Thus with the choice ((TT|l the renormalised operator iCeren = 2^kKq decouples 
from D^O. Using a standard ^-factor we find 

(KerenKeren) = -128M^ (l + gHl + ...) - traces (12) 

with the loop-counting factor M = N/ (47r^ ) . The values 

71 = 4M, 72 = -6M2 (13) 

clearly identify the operator Kq as belonging to the BMN-multiplet starting with the super- 
primary state Tr{Z^i^^) + Tr{Z^^^i) analysed in 

For minimal subtraction in the tightest sense we might choose Kg = {2, 0, 0}— 2 {1, 1, 0}. 
The divergence in Zk then sees the 0(e) term in {KqD^O) and forces a renormalisation 

Keren = Zr K - {g' ^ + + ^^o) + »(/)) D^O (14) 

where Zk still has standard form. This is not multiplicative as Zk = 1 + 2g'^M/e + . . . . 
To the order we work there is apparently freedom in the system, the stable feature being 
the anomalous dimension of Kq. 

For the spin 3 mixing problem we return to our somewhat more elegant first scheme, 
namely to make the operators tree-orthogonal at leading and next-to-leading order. This is 
achieved by 

D^,0 = {3, 0,0} +3 ({2, 1,0} + {1,2,0}) +2 {1,1,1}, (15) 
D,Kg = {3,0,0}-(l + e)({2,l,0} + {l,2,0})-(2 + e){l,l,l}, 

= (2-^){3,0,0}-(9-^)({2,l,0} + {l,2,0}) + (24 + ^){l,l,l}, 

V2 = {2, 1,0} -{1,2,0}. 

^The convention used in '25' differs by a rescaling of the coupling constant. 
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The operator V2 is the only one that is odd under reversal of the three sites in the trace 
whereby it completely decouples at classical and quantum level. The first two operators 
D^O and DzKq are defined by naive differentiation of their spin 2 ancestors. As already 
mentioned, D^O is a null-vector of the loop-level mixing; in conclusion this operator does 
not mix either. 

One may wonder what correlation functions D^Kq and Vi have w.r.t. the conformal 
primary Kq itself. By direct calculation 

{K,renD,K,ren) = -1792M^(1 + ^-^ M + . . .) ^^^^^^^l^^^^^^^^ - tvaces (16) 

as expected by differentiating formula H12() .^ More interesting is the question as to how Vi 
can be made orthogonal to Kq. We start with the definition 

Viren = ZyVi - BD.Kq (17) 

and impose 

{Vlren K^ren) = 0. (18) 

Here K^ren must be as before, while we take Zy, B in standard form ((7|, Q but with 
unknown coefficients. The simultaneous expansion of (|18j) in e and yields three conditions 
from the vanishing of the constant term at g^ and the pole and the simple logarithm at 5^. 
The result is 

Zv = I + g^-— + 0{g^), B = g^ ^ + gU—— + b2Q) + 0{g^) . 19 
4e 7 V 56e / 

In the spin 3/spin 3 mixing Viren stays orthogonal to D^Kq. The requirement of finiteness 
of {Viren ^iren) allows oue to fix the g'^ part of Zy which is in fact of the form 



Last, V2 may be renormalised multiplicatively as it does not mix. Our main result is 
that Vi, V2 have identical anomalous dimensions 



7i = ^M, 72 = -^M^ (20) 
2 16 



i.e. the values predicted in [T 



3 The graph calculation with J\f = 2 superfields 
3.1 A shopping spree in AT = 2 land 

The sl{2) sector contains the operators defined in though in general with k insertions 
of elementary fields that may all carry covariant derivatives. Independently of the spin, the 
operators are in the same [0,A;,0] representations of 5C/(4) as the 1/2 BPS states without 
derivatives. The definition is a highest weight projection in an A/" = 1 formulation, while 

^Likewise, the direct evaluation of {DzKsren D^K^ren) yields the expected a;-dependence with the correct 
normalisation 28672M^(1 + 29/14 (;^M + . . .) . The point may be worth mentioning as the projection method 
in the graph calculation is relatively involved. 
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in 7\A = 2 we may realise representatives as chains of k copies of the hypermultiplet . One 
returns to the physical scalars of the = 4 multiplet by putting 6 = 0. The TV = 2 
two-point functions of composites of matter fields involve relatively few supergraphs. The 
advantage is quite clear: at 0{g^) we have only one contributing diagram and at g'^ there 
are no more than fifteen singular graphs, which can be evaluated by standard means since 
they lead to two-point integrals up to and including three-loop level in momentum space. 

Quantum calculations for the BPS operators themselves give zero, but divergences do 
arise when the hypermultiplets (or Z's) carry Yang-Mills covariant derivatives. We find 
that the partial derivative in 2?^ = dfj_ + ig Afj_ does not introduce singularities but that 
the connection part gives rise to a set of singular diagrams. This seems logical since the 
hypermultiplet itself and its partial derivatives are short objects from an = 2 perspective 
while the connection is not. 

We will not review the N = 2 formalism in this letter; the interested reader is referred 
to jni- Our conventions are taken from the second reference therein. The N = 2 x-space 
supergraph technique was developed in the series of papers j^. The = 4 action is 

SN=i = Sq + Sv + ghost -I- gauge-fixing (21) 
Sg = -Tr{j dudC^{q+ D++ q+ + giq+[V^+,q+])) (22) 

292 A, „ V (12)(23),,,(nl) ' ^ ' 

with = {xa,G^ ,0^} and z = {x,9,9}. The "harmonics" u € SU{2) are auxiliary bosonic 
variables. We use the shorthand (23) = U2^u^-, where 2,3 label distinct points and i € 
{1,2}, while the row index is split into {+,— } according to the charge under the U{1) 
subgroup. The u-variable is used as a projector: 0^ = ufO^ and similarly for 9 and the 
spinorial derivatives D,D. The "analytic x" 

= - i(0+cj^^- + e-a''e+) (24) 

is annihilated by , Z)"*" . In the YM action the u- variables distinguish the various In- 
fields. The SU (2) structure of the three-vertex is in fact completely antisymmetric so that 
the simple product from (|23|) becomes a commutator. 

The two unconstrained quantum fields of the formalism are the hypermultiplet q^{xA, , 
0~^,u) and the YM prepotential V'^^{xa,0~^ ,0~^ ,u). We have the propagators 

ni2 = (g+(l)g+(2)) = e^^- \ , (25) 

A12 = ^ [{1-2)9+ d^et + {2-i)etdiet + e+diO^ + e+dief] , 

and in Feynman gauge 



(y++(l)y++(2)) = 6iuuu,) ^^^ %^yi (26) 



With the notation we mean 9 at point 2, projected by u'^ from point 1. (The expression 
for the propagator remains reasonable due to the presence of the SU{2) delta-function.) 
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Note that there is no higher matter /YM interaction than the cubic term in Sg, while there 
exist YM self-interactions of any order. Our set of graphs involves only the cubic YM vertex, 
though. Prom the interaction Lagrangian we retain 



J dudC^Tr{q+[V++,q-^]) , (27) 

9i (-12 , , , Tr{V++{z,u^)[V++{z,u^),V++{z,u^)]) 

— Id z auT auo du^ — — — — — : -. (28) 

12 7 I 2 i (12)(23)(31) ^ ^ 

The SSDR scheme means to leave the spinor and SU (2) algebra untouched but to modify 
the part of the propagators as in ©. 

Finally, we need to discuss the connection in = + giA^. In the A-frame ^HI we 
have 

Vt = Dt, Vt = Dt. (29) 

(Hence the + projected spinor derivatives are YM covariant.) The - projected spinor deriva- 
tives inherit a connection from the SU{2) derivative T) : 

V— = D— + giV~- (30) 
V- = [V-,D+] = D- - gz{D+V-) (31) 
■DT = [V-~,D+] = DT - gi{D+V-~) (32) 

The super-algebra does not change under covariantisation: 

{V+,VT} = -2iVo,a (33) 

from where we may read off 

V^ = + ^{D+a^D+ V—) . (34) 
The field V has a series expansion in terms of V'^^ quite like the YM action: 



V++{Z^U2) 

(12)2 



V-{z,u,) = I dn, ; + (35) 



g_^ f [y^^(.,n.),y^^(.,n3)] 
2 J ^ (12)(23)(31) ^ ^ 

We draw the reader's attention to the non-linearity of the connection (|34|) in V^^. For our 
purposes only the first two terms of the expansion are needed. Note that the linear term 
has the following two-point function with V'^^: 

( >+a,Z)+ V-m V^H2) )o,A^o = " ^'-'^'^'^'^ , ,f "^'^^ (37) 
4 ^' ^ 2tt'^[xi — XA2) 

When evaluating Feynman graphs with connection ends it is a very useful fact that the 
exponential shift in this formula coincides with one of the terms in the shift acting on 
the matter propagator. This exponential is usually reproduced also in graphs involving the 
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quadratic part in the connection because cross terms from the denominators of the attached 
propagators are ehminated when putting the outer 0's to zero. 

Let us comment on how to evaluate supergraphs. A YM-hne between two matter /YM 
vertices is relatively easy to deal with: owing to the presence of the u delta-function and the 
nilpotent numerator we may change the denominator of the V-propagator from xi — X2 to 
XAi — XA2- One SU{2) integration can trivially be done using the delta-function. We may 
now use the numerator of the propagator as a delta-function for the -f- projection of the 
two 0's, the minus projection is absent. Thus we replace, say the spinors in the exponential 
shifts A at the second vertex by those from the first and integrate away the delta-function. 
The second 9 and u integration is left over for the moment. 

Graphs involving the cubic YM vertex or a connection end are trickier in that a given 
6 may occur in contractions w.r.t. several different n's. If the YM lines end on matter 
lines it can be an advantage to shift to a chiral basis at the complicated end, otherwise it 
is hard to give general rules. In order to integrate out the 6 at the cubic vertex it is always 
advisable to choose that projector which occurs most often and to use the "cyclic identity"^ 
to write all occurrences of the 9 in question in terms of the -|-, — projections w.r.t. this one 
harmonic. 

Eventually one has to expand the exponentials and the numerators in the remaining 
0's, leading to a linear combination of differential operators under the x-integrations with 
harmonic (i.e. SU{2)) integrals as coefficients. These li-integrations can always be done 
in an algebraic fashion: if there is no denominator one can straightforwardly express them 
in terms of SU{2) epsilon tensors, else one pulls a D^~^ off the numerator and partially 
integrates it onto the denominator of the integral thus producing delta- functions UH]- For 
the more complex g"^ supergraphs we admittedly had to resort to a Mathematica programme 
to do parts of the Grassmann algebra and the SU{2) integrations. 

In the pictures below we have split the lines emanating from the outer points. We do 
so purely for convenience of drawing — all lines coming in from the left start at zi = 
{xi,9i,ui} and all lines ending on the right come together at the point with coordinates 
Z2 = {x2,92,U2}- "Spectator lines" , i.e. free matter propagators have been omitted. There 
are, of course, as many as needed to connect all matter fields in the composite operators 
at points zi, Z2. We did not indicate the outer partial derivatives either since they do not 
interfere with the Grassmann- and 5C/(2)-integrations. 

Note that any two-point graph related to a [0, k, 0] operator will come with a factor (12)*^, 
which we simply drop. It is equal in the tree-level and the and g'^ parts and universal 
to all operators built from k fields; therefore it does not infiuence the renormalisation at 
all. At order g"^ individual graphs also produce "harmonic-nonanalytic" terms coming with 
(12)^+^(1^2"), which must eventually cancel. A good sign is that all such terms contain 
free propagators times the finite integral / from graph Bq below. We did not explicitly check 
the vanishing of their sum because they are irrelevant to the second anomalous dimension. 

^For three harmonics: (12)w3 + (23)mi + {31)u2 = 0, shnilarly with one or several - projections. 
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3.2 Graphs 

At order there are only two graphs, for which we find the x-space expressions 



So = -(12)2 □!/(!, 2; 1,2), (38) 
Go = -(12) idi^g{l,l,2). (39) 

Here we have omitted the normalisation constants of the propagators, i.e. cq = — l/(47r2) 
from the hypermultiplet propagator and — 2co from the normalisation of the y-line. The 
integrals /, g are as the pictures suggest — there is no extra numerator and all lines are to 
be replaced by l/a;?^-, hence 

^' = / .2 t^"f"J^2 > 9{l. 1,2) = f . (40) 

The topology of a graph cannot change due to the spinor and SU (2) algebra other than by 
going into a "derived topology" , i.e. one in which some lines are shrunk to points by box 
operators. 




Figure 1. Graphs at order g"^ 



Let us focus on diagram Bq. The double-box integral / is finite and by dimensionality 

we have /(I, 2; 1,2) oc l/x\2i so that the outer box operator prodviccs a delta-function. 
In the SSDR scheme this is seen as an 0(e) term, see Section 2. This would be the only 
contribution in a two-point function of products of hypermultiplets without any derivatives. 
Indeed, such operators are 1/2 BPS and their two-point functions should not receive any 
quantum corrections except for contact terms. 

What happens when the diagram has partial derivatives on the outer legs? The Min- 
cer system can calculate / in momentum space, and it can also handle arbitrary scalar 
numerators. We use the programme according to the following strategy: 

• All partial derivatives are translated into momentum space as momenta with open 
indices. We make the sets of indices symmetric and traceless separately at the left 
and at the right end of the graph, so that they define spin si and spin Sr "harmonic 
tensors" , respectively. 

• There are Min{si, Sr) + 1 independent ways of projecting this with the total ingoing 
momentum q, because we may choose to contract through Min{si,Sr) indices with 
the metric ij. All remaining open indices are contracted on q. 

• For any given distribution of outer partial derivatives/momenta we calculate all such 

projections by Mincer. 

• The integral with open indices (symmetrised and traceless at each end) contains ex- 
actly as many independent Lorentz structures as there are projections, the criterion 
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being once again how many r/ symbols with one index from the left set and one from 
the right set are involved. It is a simple exercise to set up and invert the set of linear 
equations which allows to reconstruct the coefficient of each basis element from the 
projections. The coefficients are generally e-dependent; we keep 0(1), 0(e). 

• We Fourier-transform back to j;-space. 

In this context, too, it is actually most convenient to discard all traces since under the FT 

^ ''^'-'J^-^ j^j^ (41) 

and the term without rj is simply obtained by having all derivatives act on the denominator. 

For Bq we have scanned all spins up to three derivatives at one end and four at the other. 
Our findings are surprising: the term without 77 remains a contact term for any harmonic 
tensors at the left and at the right end, while terms with ij may be finite. In conclusion, 
the "BPS-like" diagram Bq drops out of the calculation of the first anomalous dimension. 

The calculation of Go can be done in a variety of ways because the underlying integral 
g{l,l,2) is elementary. Since our computer programmes were originally geared at the [0,2,0] 
operators we have chosen to join one free line and compute the resulting two-loop diagram 
as a derived topology of the T2 graph f2T|, alternatively one may use Tl. The results of 
Mincer have to be scaled down by a factor (47r)^' where I is the loop order of the given 
integral. In configuration space graph Go has a simple pole.^ 

We remark that the diagrams have to be recomputed for every different distribution of 
derivatives. The combinatorics becomes more elaborate, too, since equal looking legs can 
be distinguished by the number of partial derivatives. 

On the computing side, we have implemented the large combinatorics for the [0, 3, 0] 
operators in a Mathematica programme, which uses a table of integrals to calculate the 
matrices of two point functions at tree- level and orders g'^,g'^- The calculation of each inte- 
gral requires the construction of the necessary harmonic tensors and their projections; for 
convenience of programming we still do this under Mathematica. Finally, we link in Mincer 
under Form to evaluate the projected integrals. 

Let us switch to order g^. There are five "BPS-like" graphs, i.e. topologies that one also 
finds in two-point functions of 1/2 BPS operators, though without extra partial derivatives 
on the outer legs. When analysing two-point functions of 1/2 BPS operators one tends 
to find two linear combinations of graphs coming with independent colour structures j24j . 
Protectedness is tantamount to both of these linear combinations being contact terms. In 
the manifestly finite setup we use {M = 4 in terms of = 2 fields) the only divergences are 
associated to composite operators. Therefore, graphs related to BPS objects should even 
be individually finite. 

The BPS-like supergraphs are drawn in Figure 2. In B2 there is no vertex in the middle 
of the diagram, likewise in B^ the YM line between points 5 and 6 does not touch the matter 
line in the middle of the diagram. 

®The index of the connection "migrates" to the opposite end of the graph, see (1391 . When building up 
the harmonic tensors we must put this index into the tensor at the left end, of course. 
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Figure 2. BPS-like graphs at order g'^. 



We find 



Bi = {UfniDi LA, (42) 

B2 = (12)2 [-□i(9i3 + 524)' + (i2)(i-2-)(ni5n24 + ni3n26)]iv^', (43) 
^3 = (12)2 [-□in34 + (i2)(i-2-)(ni5n23 + ni3n26)] ^^^^ 

54 = (12)3 [(5i3 + 914)2(525+526)' -□in45/2-ni6n23+ (45) 

+ (i2)(i-2-)ni6n23] ^4, 

where the normahsation constants of the propagators were again omitted. The x-space inte- 
grals LA, NP, BE are the so-called "ladder" , "non-planar" and "Benz" three-loop integrals 
as depicted in Figure 2. We illustrated by equation (|4()|) how to read them off from the 
graphics; there is no additional numerator and the denominator is given by the appropriate 
product of l/xfj, no propagator is doubled. By dij, Oij we mean derivatives acting at point 
z on a line going from i to j, whereas Di means an outer box-operator acting on the whole 
expression. 

All three topologies LA, NP and BE are finite without a numerator, but differentiation 
under the integral could potentially make them divergent. Using the Mincer system we have 
checked that the complete expressions Bi, B2, -B3 do not diverge either. Like Bq the graphs 
remain finite, too, when symmetrised, traceless products of derivatives are introduced acting 
from the left and from the right, respectively. A scan found no singularities up to spin 3 at 
one end and spin 4 at the other. 

The integral in graph i?4 was baptised "P4" above as it is planar and has four loops. A 
quick evaluation of ()45() in point-splitting suggests that the graph without extra derivatives 
is finite. Unfortunately, we have no way of calculating B4 in dimensional regularisation, in 
particular in the presence of outer partial derivatives. By analogy with Bq . . . B^ we assume 
that B4 is finite, too, so that we may dispose of it. Diagram B^ should behave accordingly; 
at any rate, as a non-planar structure it is not relevant in the context of our leading 
analysis. 

Note that the non-analytic parts in (|43() ... (|45() always reduce to /(I, 2; 1,2) whence 
they are finite by what was said about Bq. 
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The singular graphs at order g'^ arc hsted in Figure 3. We have omitted the mirror 
graphs and some diagrams that must carry at least one power of rj on dimensional grounds, 
e.g. the graph with three connection lines joining at a cubic YM vertex. The product 
Bq * Go has no singular term without rj either, see above. 

The graphs in the third row carry an index /v, at point 1 on the left and an index u at 
point 2 on the right. The dot in G4, G-j, Gg indicates the quadratic term in the connection 
from which two F-lines emanate, whereas G13, G14, G15 have two different occurrences of 
the connection on the left with the indices fj,,!^. 




Figure 3. 0{g'^) graphs with connection lines. 



Evaluation of the Grassmann- and S'L''(2)-integrations reduces the supergraphs to the 
following differential operators under the appropriate integrals: 

Gi = {12fi[ 8^,1034/2 - d^s^Di] (46) 

G2 = (12)2i[ 9^in35/2 + 5^23(^14 +5i5)'] (47) 
G3 = il2fi[ 5^in35/2 + 5^23(5l4 + 9l5)'-5^23ni4/2 (48) 

- 5m1°45/2 - 5^24 (5l3 + 5i5)2 + a/,24ni3/2 

- (5^23 - 5^24 )□ 15/2 

- {df,23 - 5^24)ni/2 

- 5^i3(ni + 023 - 2n24)/2 + d^uiOi - 2^23 + □24)72 

+ 5^1 (023 - ^24) - (12)(l-2-)9^i5(n23 " ^24)] 
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G4 


= (12)^^[ 


-(5^23-5m24)/2] 


(49) 


Gs 


= (12) ^[ 


-a^2ni5-(i2)(i-2-)a^i5ni] 


(50) 


Ge 


= (12) i[ 


- 9^2ni5 + 5^2ni3/2 - d^a^i + (12)(l-2-)5^i5ni] 


(51) 


Gj 


= (12) i[ 


9,2] 


(52) 


Gs 


= (12) i[ 


9,2] 


(53) 


Gg 


= (12) [ 


- d^34du34: + rj^^n^i/A] 


(54) 


Gio 


= (12) [ 




(55) 






-{l + {U){l'2~))d^ndu2A] 




Gn 


= (12) [ 


— dyi3d^2'i + dnidu2 — (9^14^^1/23 + dfj_i3dy2A)/'^ 


(56) 






+ r/^.(9l3 - 523). (5l4 - 924)/4 - (12)(l-2-)a^i4a,24] 




Gl2 


= (12)M 


d,2?,dviA ] 


(57) 


Gi3 


= (12) [ 


d^34du2 + ^,1^034/4.] 


(58) 


Gi4 


= (12) [ 


d^23dlu - du23d'/ii4 - {du23d'f,i4 - dfj,23d'li4)/2 


(59) 






+ V,uidl3-d23).{d[4-d'{4)/A] 




Gi5 


= (12)^ [ 


— 9^23^^1/24] 


(60) 



The combinatorics for G4 is like for G3 because the quadratic term in the connection is 
a commutator. The graphs come in a fixed sum in which G4 compensates the Dis term 
in the expression for G3. The box operators □23,024 in (|48|) reduce the underlying BU 
topology to the usual finite /(I, 2; 1,2). In all these terms the remaining derivative acts 
from the left and thus there are no poles. The BU integral itself is in fact finite, too, as 
long as no double lines are introduced. We have checked that the Di terms in the fifth line 
of H48() may be dropped as well. The integral is slightly more robust than / in that the 
term (9^23 — <9^24)ni is not singular either, although the derivatives act from the right. In 
conclusion, 

G3 + G4 = {I2fi[ d^iD35/2 + d^23{du + di5f -d^23Du/2 (61) 

- 045/2 - a^24(9i3 + di^f + a^24ni3/2 + . . . ] 

where the dots indicate the omitted finite terms. Similarly, G5, Gg, Gj come in a fixed sum. 
On discarding finite terms: 

G^ + Ge + Gj = (12) i [ - d^2Di5 + 9,.2ni3/2 + • . .] (62) 

(The numbering of the points refers to G^,Gq.) Owing to the box operators both terms 
can be fitted into the 02 topology after joining a free line. It is interesting to note that the 
□13, 014 terms in H61() are of the same type, while the remaining terms are as in G2. The 
graphs Gi . . . G7 seemingly always sum to zero if one end of the correlator is a BPS state. 
Diagram Gg has the same combinatorics as the other graphs in the second line, but it is a 
separate piece in this sense. 

In diagrams Gio,Gii the potentially divergent pieces without rj are 

Gio = (12) [ 9,,149^23 + 9^349i.34 + •••], (63) 

Gn = (12) [-a,i3a^23 + ...] ■ (64) 
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In our calculations the complete sum of graphs turned out to be finite within each class. 

One last simplification concerns graph G14: for our purposes we may once again drop 
the T] term. Further, when the two indices are symmetrised 

Gi4 = (12) [ 9(^23 {Kn, -d'l^,,)/2 + ...]. (65) 

We remark that diagram G14 starts to contribute only at spin 4. 

All the integrals can be evaluated by the Mincer system just as the order g'^ examples: 
where there is only one matter line we join a second, free one. In this way Gi, G5 . . . Gg, 
G12 ■ ■ ■ Gi4 are of topology 02 (or derived of 02) while the remaining diagrams belong to 
the BU category. The rest of the algorithm is as before, although formula (|41j) must be 
adapted to the dimension of the three-loop integrals: 

(.,)(--)(,^)-- ^ ''^'-'J^---^ (-ai)(---) -^-^ (66) 



4 Conclusions and outlook 

We have discussed the descendant structure and renormalisation of the sl{2) sector operators 
of length three up to spin three. These are a 1/2 BPS state O, a spin two operator Kq 
belonging to the second lowest two-impurity BMN multiplet and two spin three primaries 
Vi, V2 with identical scaling dimensions up to order g'^. The anomalous dimensions are 

15 225 

ife : 71 = 4M, 72 = -6M^ , Vi, V2 : -fi = — M, 72 = -—M^ 

2 16 

with M = N/{4iT^). The result for Vi, V2 is in exact agreement with c.f. Table 2. 

The most difficult step in developing any higher-loop dilatation operator is to make it 
work when the spin-chain length is smaller than the maximum number of adjacent sites on 
which the operator can act ("wrapping"). From this perspective, the hardest test on any 
universal formula is to go from the shortest spin-chain length to operators with subsequently 
more sites in the trace. In the sl{2) sector this would be the step from twist two to twist 
three, which we have checked for the Bethe-ansatz 16„, initially at the lowest non-trivial 
spin. 

The obvious continuation of our work is to investigate other mixing examples at higher 
spin as well as for operators with more elementary fields. The ultimate aim of the line of re- 
search is to set up a two-loop dilatation operator for the sl{2) sector. It would be interesting 
to see whether the dilatation operator can be directly constructed from the elements of our 
renormalisation procedure by going back to the non-orthogonal basis {si, S2, S3} defined in 
equation (P). 

The renormalisation of Vi is non-multiplicative: we have put Viren = ^1 ~ 13 DzKq 
with a singular factor B that cannot be split into h{g'^)Zy, where h{g'^) is some finite 
function. This example of orthogonalisation in dimensional regularisation closely parallels 
the discussion of the Konishi anomaly in j25j . 

Last, the limiting step in our current software is the manipulation of harmonic tensors, 
because for higher spin these involve very many terms. It would clearly be desirable to 
transfer a larger part of the calculation from Mathematica to Form. 
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